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The WMAP haze is an excess in microwave emission coming from the center of the Milky Way 
galaxy. In the case of synchrotron emission models of the haze, we present tests for the source 
of radiating high-energy electrons/positrons. We explore several models in the case of a pulsar 
population or dark matter annihilation as the source. These morphological signatures of these 
models are small behind the WMAP Galactic mask, but are testable and constrain the source 
models. We show that detailed measurements of the morphology may distinguish between the 
pulsar and dark matter interpretations as well as differentiate among different pulsar models and 
dark matter profile models individually. Specifically, we find that a zero central density Galactic 
pulsar population model is in tension with the observed WMAP haze. The Planck Observatory's 
greater sensitivity and expected smaller Galactic mask should potentially provide a robust signature 
of the WMAP haze as either a pulsar population or the dark matter. 

PACS numbers: 95.35. +d, 98.70.Vc 



I. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
has provided a detailed map of the cosmic microwave 
background (CMB) as well as the foreground from our 
Galaxy [TJ [2] . The foreground emission from the center 
of the Galaxy shows an excess of emission in the inner 5 — 
20° around the Galactic center. This excess of microwave 
emission is known as the "WMAP haze" [3]. Due to 
the Galactic plane, the current data on the haze only 
goes down to about 6 degrees from the Galactic center. 
Two notable features of the haze are the approximate 
spherical symmetry and the strong angular dependence of 
the flux. In particular, the flux increases quickly toward 
the Galactic center. 

The size and shape of the haze is currently not well con- 
strained. Using a different foreground model may change 
the look of the haze slightly. For instance, Ref. jl] used 
a different synchrotron map to model the WMAP data 
and found that the haze in that case is smaller in spa- 
tial extent and slightly smaller in magnitude. A cleaner 
foreground subtraction is needed to further constrain the 
source of the haze. In this paper, we will use the mea- 
surement of the haze presented in Ref. [S]. 

An explanation for the WMAP haze is a previously 
unknown source of microwave emission. The frequency 
dependence of the haze is quite hard, so a hard source 
like synchrotron radiation is a likely candidate [5H7] . The 
magnetic field in the Galactic center tends to be tens of 
microgauss. Therefore, the synchrotron radiation from 
highly relativistic electrons and positrons near the Galac- 
tic center could be the source of the microwave signal [H]. 
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Another possibility for the source of the haze was ther- 
mal bremsstrahlung from ionized gas [3j. However, the 
Ha skymap shows no significant increase in the regions 
where the haze is strongest. At high density and high 
temperature (~ lO^K) such a gas could still explain the 
haze, but the emission from such regions is constrained 
to be small, ruling it out as the source of the haze [S] • 

In the synchrotron emission interpretation, calculating 
the complete propagation of electrons in the interstellar 
medium requires the full diffusion-loss equation. This in- 
cludes spatial diffusion of the e+e~, reacceleration of the 
particles due to momentum-space diffusion, energy loss 
due to several different mechanisms, convection of the 
particles in the Galaxy, and the particle source. To solve 
this complete equation, one may use the software GAL- 
PROP [inilTT]. However, the major necessary physical 
features of the source and synchrotron emission can be 
sufficiently modeled using a Green's function solution to 
the diffusion- loss equation, which we employ [T^ ITB] . 

The source of high energy electrons and positrons re- 
quired for synchrotron emission remains a mystery |14j . 
One intriguing explanation of the haze is from dark 
matter annihilations in the Galactic center [B]. The 
e+e^ produced from by the annihilation products move 
through the Galactic magnetic field, creating the diffuse 
synchrotron emission [S]. This model matches the haze 
well, especially the strong angular dependence and ap- 
proximate spherical symmetry. 

It should be noted that Ref. 15J claims that dark 
matter annihilations cannot explain the haze because it 
would require too large of a dumpiness boost factor. In 
this work, we employ similar methods to those used in 
that paper and do not find that too large of a boost factor 
is necessary. Ref. |15j averaged the dark matter density 
over the direction perpendicular to the Galactic plane, 
which is inaccurate in light of the morphological effects 
we present below. 
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Another source for the haze is e+e~ coming from the 
magnetosphere boundary of Galactic pulsars [16]. This 
type of emission from nearby pulsars can be the source 
of positrons in the PAMELA results [T7H1I] and also the 
source of features seen by ATIC and the Fermi 

Gamma-Ray Space Telescope [23j [24] . 

In this paper, we explore models of the WMAP haze 
in the context of synchrotron radiation from the elec- 
tron and positron production of Galactic pulsars and 
dark matter annihilations. We will show tests of this 
diffuse synchrotron emission which could distinguish be- 
tween the exotic explanation of dark matter annihila- 
tions and the astrophysical pulsar sources. We also show 
how haze observations can test pulsar and dark matter 
models separately. These tests can be done with upcom- 
ing results from the Planck Observatory, which has been 
shown, e.g. in Leach, et al. [25 , to be expected to have 
a much smaller Galactic mask than WMAP [26l [27] . 

Recently, Kaplinghat et al. [16^ showed that the mor- 
phology of the haze is different depending on the source. 
This work pointed out that for a source with spher- 
ical symmetry, the haze should be slightly elliptically 
stretched along the Galactic plane while for a centrally- 
peaked pulsar source, the signal lies primarily along the 
Galactic plane and has less spherical symmetry. 

In this work, we examine the morphological structure 
of the WMAP haze in detail. First, we consider the signal 
due to pulsars, both from a Gaussian distribution and one 
which vanishes at the Galactic center. Second, we look at 
dark matter annihilations models' sensitivity to the dark 
matter density profile. Lastly, we compare the pulsar and 
dark matter scenarios to quantitatively distinguish them 
from one another. 



II. CALCULATION OF SYNCHROTRON FLUX 

To calculate the propagation of the electrons from 
their source, we must solve the diffusion-loss equation 
for e+e^: [T1[TH] 
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diffusion-loss equation is derived from the full equation 
by assuming that reacceleration is small compared to spa- 
tial diffusion and that convection is small compared to 
the particle velocities. For highly relativistic e^e" in the 
Galaxy, these are both reasonable assumptions. 

For the simplified diffusion- loss equation, we shall use 
the parameterization of Ref . [5] . We assume a spatially- 
independent diffusion constant K [e) = lO^^gCi-SS cm^s~^ 
inside a "diffusion zone" of half-thickness 3 kpc away 
from the Galactic plane and "large" in the radial direc- 
tion of the Galactic plane. Outside this region, it is as- 
sumed that all particles are no longer confined to the 
Galactic magnetic field and therefore are free-streaming. 
Inside this region, the Galactic field is assumed to be a 
constant. This is a reasonable approximation because 
it is known that the Galactic magnetic field decreases 
quickly away from the Galactic center, but its exact 
profile is unknown. The energy loss rate is also as- 
sumed to be spatially-independent and is parameterized 
by 6(e) — /te with r^; = 2 x 10^^ s. This characteris- 
tic energy-loss time te accounts for losses due to inverse 
Compton scattering on the CMB and starlight, as well 
as synchrotron energy losses. The source term Q {E, x) 
is the rate of particles created per unit time per unit vol- 
ume per unit energy. We will further assume that the 
diffusion- loss equation is in steady-state, so all terms are 
time- independent. The assumptions of spatial and tem- 
poral independence of the diffusion and loss parameters 
should be valid as long as we are concerned only with re- 
gions close to the Galactic center, so the average values of 
parameters can change only slightly in the time and dis- 
tances involved. The assumption of time-independence 
in the source term assumes that the rate of pulsar forma- 
tion and death is steady over the timescale of the electron 
propagation, during which only a few new pulsars should 
have formed. 

For our parameters, an analytic solution to the 
diffusion- loss equation can be found using Green's func- 
tion techniques [TH |T3]. Writing the source term as 
Q{E,x) = f{E)g{x), this analytic solution is given by: 



dn 
dE 



nOC 

E-^ / dE'f{E)TD{E',E,x) , 

J E 



(2.2) 



with K the diffusion constant, h the energy loss rate, Q 
the source term, and e = E / {\ GeV). This simplified 
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where x = {s, 0, z) are cylindrical coordinates from the 
Galactic center, i = 3 kpc is the end of the diffusion zone 
and D'^ = 8.4 (e-"-67 _ g'-o.67^ jg ^j^g integrated 

energy-dependent diffusion parameter. To convert this 
e^e~ density to a synchrotron flux at earth, we must 
integrate over line-of-sight ^ and then convolve it with 
the synchrotron spectrum for a single electron |12j : 
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Here, Tgyn = 4Tb is the energy loss time due to syn- 
chrotron emission, y = {vm^c^) / {vbE"^), v is frequency, 
and vb — (ieH) j {^'nrrig) is the characteristic synchrotron 
frequency. As in Ref. [S], we choose a constant magnetic 
field of 10 /iG in the Galactic center and look primarily 
at the V — 22 GHz band of the haze. We will use a coor- 
dinate system where a is the angle out from the Galactic 
center and t\} is the angle up from the Galactic plane. 



III. ELECTRON SOURCE FUNCTION 

A. Source Function for Pulsars 

In order to calculate the diffuse synchrotron flux <!> 
coming from pulsars, it is necessary to solve for the rate 
of e+e~ produced per units energy, volume, and time. 
This can be separated into two parts, the rate of electron 
(and positron) production in a given pulsar, d^N/ {dEdt), 
and the spatial distribution of the pulsars in the Galaxy 
p{x). 

For the number of electrons and positrons produced at 
a given energy in a given time, we use a model which 
does a fit to gamma-ray production from Galactic pul- 
sars to determine the e^e^ rate from pulsars in the 
Galaxy [28]. Above about 1.5 GeV, this model is sim- 
ilar to the positron model in Ref. [23] . For the combined 
e^e^ production rate, this model gives 
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where /e is the ratio of e+e^ to photon production in 
mature pulsars and A^ioo is the rate of pulsar formation 
in the Galaxy per hundred years. The free parameters 
/e and -/Vioo are to account for the fact that the fit to 
gamma-ray pulsar data was done assuming a pulsar pro- 
duction rate of one per hundred years, so the overall nor- 
malization is unknown. 

The distribution of pulsars in the Galaxy can be broken 
into two parts, a radial component ps (s) and a compo- 
nent against the Galactic plane pz {z). The full pulsar 



density is, then, p{x) — ps{s) Pz{z). The perpendicu- 
lar distribution pz should be peaked about the Galac- 
tic plane to match pulsar surveys with pz {—z) — pz {z). 
This has been parameterized as [3D] 
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with the normalization dz — 1. 
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Because of observational bias in pulsar surveys, there 
is much model dependence in the shape of the radial dis- 
tribution of pulsars ps- For example, one must be careful 
in selection effects due to bright pulsars being more easily 
observable than dim ones. Also, in the inner part of the 
Galaxy, the pulsar distribution is largely unknown due to 
scattering and obscuration, so a targeted high-frequency 
radio survey would be needed to detect these innermost 
pulsars [3T] . 

This ambiguity in models has lead to two distinct fits to 
the pulsar survey data: a Gaussian fit with a peak at the 
Galactic center or a Boltzmann-like distribution which 
goes to zero at the Galactic center. As evidence for the 
latter, the Effelsberg 5 GHz Galactic center survey has 
not found any pulsars towards the Galactic center, but 
it is possible that any survey below 10 GHz cannot avoid 
scattering and obscuration effects [32135]. If a large num- 
ber of pulsars in the inner kiloparsecs are observed by the 
Fermi Gamma-Ray Space Telescope, it could distinguish 
between these two models as well [23] . 

The zero central density distribution is one in which 
the function peaks toward the Galactic center and then 
goes to zero at the center itself. This model is plausible 
if the pulsar distribution is similar to the distribution of 
supernova remnants and some astrophysical gas popula- 
tions, which are in rings peaked about 4 kpc from the 
Galactic center [31] I3m36j . This model of the pulsar 
radial distribution can be parameterized as [32] 
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where s is the distance from the Galactic center in the 
Galactic plane (in kpc). This distribution has been nor- 
malized to unity using 27r ps{s)sds = 1. 

Another radial distribution of Galactic pulsars is Gaus- 
sian, with a peak at the Galactic center falling off in the 
outer Galaxy. This model is motivated by stellar evo- 
lution arguments where the distribution is expected to 
be more densely peaked near the Galactic center, follow- 
ing the stellar distribution in the Galaxy. Such a radial 
distribution has been parameterized as |30j 
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where s is the distance from the Galactic center in the 
Galactic plane. This distribution has been normalized to 
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FIG. 1. Angular flux profile $ versus angle from the Galactic center a. The diagram on the left is for the Gaussian pulsar 
distribution and the diagram on the right is for the pulsar distribution with no pulsars at the Galactic center. From top to 
bottom, the lines are the 1/1 = 0° (dashed), il> = 30° (dash-triple dotted), -i/; = 60° (dotted), and ip = 0° (solid) plot. The plots 
are scaled with the = 90° plot normalized to the flux of the WMAP haze. The points are the WMAP haze with errors from 
Ref. [5]. 



unity as the zero central density distribution was. An ex- 
ponential fit like the one used by Ref. \W can be modified 
to be indistinguishable from the Gaussian form (30j . 

B. Source Function for Dark Matter Annihilations 

To compare the pulsar source scenario with dark mat- 
ter ones, we consider two different Navarro- Frenk- White 
(NFW) [37] type dark matter profiles with WIMP anni- 
hilations. In this model, the relativistic e+e" are formed 
when dark matter particles annihilate in the Galactic 
dark matter halo. This can happen if the dark matter 
is its own antiparticle, as in the case of supersymmetric 
neutralinos [5]. For the halo shape, we use a generaliza- 
tion of the standard NFW profile ^37] . The generalization 
can be written in the form 

^ " (r/r,)7(l + r/r,)3-7 ' ^^'^^ 

for scale density po ^nd scale radius rg. The origi- 
nal NFW profile has 7=1. A slightly steeper profile 
with 7 = 1.2 is also consistent with Milky- Way con- 
straints [SHHID]. We take — 25 kpc as out canoni- 
cal scale radius. The scale density is normalized to lo- 
cal earth dark matter density (~ 0.3 GeV/cm^) [41]. A 
standard WIMP has a weak-interaction cross-section of 
(crti) = 3 X 10~'^^ cm^s~^ [5JI12]. For simplicity, we con- 
sider the direct annihilation channel where the product is 
exactly one positron and one electron with energy equal 
to the dark matter particle mass. We use 100 GeV dark 
matter particles as a canonical value. The morphological 
tests here are not sensitive to this choice. The source 



function for e^e formed from such a dark matter anni- 
hilation is 

QiE,x)^'^ = ^{<JvmE M) . (3.6) 



IV. SYNCHROTRON FLUX FROM SOURCES 

A. Flux from Pulsars 

For the standard pulsar distribution which increases 
toward the Galactic center, the e+e" source term is 

Using this with an overall scaling factor for the normal- 
ization of the particle production rate we find that the an- 
gular flux profile (<& vs. a) away from the Galactic plane 
fits the WMAP haze well, in agreement with Kaplinghat 
et al. [TB]. As in that work, we use an overall scaling fac- 
tor to account for the uncertainties in particle production 
as well as a constant offset to account for uncertainties 
in background subtraction. We flnd that a physically 
reasonable scaling factor with electron production with 
similar magnitude to photon production in pulsars fits 
the data well. 

We find that if one measures at an angle not orthogonal 
to the Galactic plane, but in an arbitrary direction, there 
is a strong dependence in the angular fiux profile. As the 
direction of the line-of-sight shifts away from the vertical, 
there is a significant flattening in the angular flux profile 
(Fig. [T|). This leveling-off is due to the large scale with 
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FIG. 2. Relative flux $ versus Galactic coordinates i and h for the Gaussian pulsar profile (left) and the pulsar profile with no 
pulsars at the Galactic center (right). The mask is that from WMAP K-band mask Kp4 [43] . 



which the pulsar distribution drops off in the s-coordinate 
as compared to the relative steepness with which it drops 
off in the z-coordinate. Also, the overall magnitude of <I> 
changes by a factor of four to five at angle a = 10°. 
This effect is due to the fact that all fluxes should be 
equivalent in the Galactic center and begin dropping off 
with varying steepness at that angle. The 2-dimensional 
flux profile is shown in Fig. |2] 

For the zero central density pulsar distribution, the 
source term becomes 



dEdt 



(4.2) 



With this profile, we need a slightly different normal- 
ization to match the fiux of the haze, which could be 
accounted for by modifying diffusion- loss parameters. In 
this case, however, the fiattening-out effect happens at 
angles much closer to the normal to the Galactic plane. 
Because this profile has its maximum density at around 
4 kpc out from the Galactic center, the line-of-sight in- 
tegral keeps increasing out to a « 25°. This creates a 
second peak in the angular fiux profile at a « 20° before 
the flux begins to fall off at large angles. This can best 
be seen in the (£, &) diagram. Fig. |2] around (±20°, 0°). 
There the magnitude increases slightly; this feature is 
not apparent in the WMAP haze. Unless this feature 
is detected, it is improbable that such a Galactic pulsar 
distribution could cause the haze. This could be seen in 
upcoming Planck observations which will likely have a 
smaller mask 1251. 



B. Flux from Dark Matter Annihilations 

Here we consider Dark Matter annihilations as a source 
for the haze. Using the source term Q{E,x)^^ and nor- 
malizing the flux $ away from the Galactic plane to the 
haze, we solve the diffusion-loss equation to determine 
the angular flux profile for the 7 = 1.0 and 7 = 1.2 
NFW-type profiles. 

The results of the angular flux proflle for these cases 
are shown in Fig. [3] As with the flux proflles due to 
pulsars, we used a constant offset to account for uncer- 
tainties in background subtraction as well as a boost fac- 
tor to account for the possible dumpiness of the dark 
matter. We found that it only required reasonable boost 
factors of five to ten to match the haze. Due to the sig- 
nificant lack of dependence on direction, we have only 
plotted the t/j — 0° and ijj — 90° (<i>, a) curves for each 
case. Significantly, the turning-over effect seen with pul- 
sars is not seen here. Towards the Galactic plane, there 
is a slight flattening of the curve, but it is not as dras- 
tic as with the pulsar source. Given proper foreground 
subtraction, this distinction may be used to distinguish 
between the dark matter and pulsar explanations of the 
WMAP haze. In the diagrams of relative magnitude of 
<i> versus £ and b, Fig. |4j the two dark matter profiles have 
very similar characteristics. Both are significantly closer 
to spherical symmetry than the pulsar diagrams. There 
is some broadening in the 6-direction, due to the size of 
the diffusion region being finite away from the Galactic 
plane. The cuspier 7 = 1.2 profile loses magnitude much 
more quickly away from the Galactic center, so it appears 
slightly more spherically symmetric on a plot weighted 
by magnitude. However, both fall off in similar ellipti- 
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FIG. 3. Angular flux proflle, $ versus angle from the Galactic 
center a for the 7 = 1.0 and 7 = 1.2 dark matter profiles. 
The upper lines are for the tp = 90° (solid) and 1/' = 0° 
(dotted) curves using 7 = 1.2. The lower lines are for the 
= 90° (dash-triple dotted) and ij) = 0° (dashed) curves 
using 7 = 1.0. The plots are scaled with the i/j = 90° plot 
normalized to the flux of the WMAP haze. The points are 
the WMAP haze with errors from Ref. (5j. 

cal regions. This sensitivity of the WMAP haze to the 
dark matter density profile can be used to constrain the 
Galactic halo's profile if the dark matter model is verified. 
In particular, this should be seen in Planck observations 
using a likely smaller Galactic mask [25 . 

C. Dark Matter versus Pulsars 

Not only can the details of the pulsar or dark mat- 
ter distributions be tested separately, but these tests can 
distinguish between the two models and determine the 
source of the WMAP haze. For comparison, we will use 
the best-fitting candidate from each category: the Gaus- 
sian pulsar distribution and the 7 = 1.2 dark matter dis- 
tribution. Within the known values for each source, both 
can create a synchrotron signal large enough to cause the 
haze. Both sources have increased flux toward the Galac- 
tic center which falls off quickly away from it. 

The largest difference between these two possible 
sources is the strength of the signal towards the Galac- 
tic plane. For the pulsar source, the strongest signal 
comes along the Galactic plane and does not decrease 
very quickly near it. The dark matter source has approx- 
imately equal signal strength in all directions, though it 
does lose strength more slowly along the Galactic plane 
than in other directions. Unfortunately, the Galactic 
plane itself masks most of this primary difference, though 
there are regions near the Galactic center which remain 
unmasked. We therefore compare the signal away from 
the Galactic plane. 



As can be seen in Fig. |5] there are two significant dif- 
ferences between the pulsar source and the dark mat- 
ter source when looking above the Galactic plane. For 
i/j > 30° , there is little distinguishable difference between 
the dark matter signal in any direction, so it is simply 
represented by the average over angles above 30 degrees. 
Primarily, the pulsar source has a different magnitude at 
different angles V'- This can be seen by the vertical shift 
in the plots for different tp values. Even among the angles 
above the mask, there is a factor of 3 decrease in the flux 
$ for lines-of-sight away from the Galactic plane. Secon- 
darily, there is a flattening in the pulsar signal that is not 
present in the dark matter signal. Over the angles a = 6° 
to 15°, the tp = 90° plot changes by roughly 50 percent 
more than the ip = 30° plot. Such an effect is seen in 
the dark matter plots as well, but below the 10 percent 
level, which changes the signal minimally. Even above 
the mask, the two sources have substantial differences. 
With the Planck observatory's smaller expected Galactic 
mask, these differences should be even more apparent. 

V. CONCLUSIONS 

High-energy electrons/positrons moving in the Galac- 
tic magnetic field could be the source of the WMAP haze. 
Using the diffusion-loss equation, we have taken a set of 
models of likely sources and calculated the resulting syn- 
chrotron signal as a test for the models. The simplified 
equation can be solved analytically using a Green's func- 
tion approach, where it is useful for an understanding of 
the underlying physics. We considered tests of several 
models of the haze source, namely a peaked Gaussian 
distribution of pulsars, a distribution of pulsars peaked 
away from the Galactic center, and two dark matter dis- 
tributions as possible sources, all consistent with other 
constraints. 

We find that the WMAP haze could be caused by dif- 
fuse synchrotron emission due to pulsars, in agreement 
with Kaplinghat et al. [TB]. Moreover, we find that if the 
haze is caused by pulsars, the angular flux profile from 
the Galactic center should be peaked more sharply when 
using a line-of-sight away from the Galactic plane and is 
fiattened significantly when using a line-of-sight closer to 
the Galactic plane. Also, if the pulsar distribution van- 
ishes in the Galactic center, a brightness about 20° out 
along the Galactic plane should be visible. 

With annihilating dark matter, the angular flux proflle 
from Galactic center is largely independent of direction 
and has rough spherical symmetry. If the signal is found 
to be strongly spherically symmetric, then this would be 
an indication that annihilating dark matter could be the 
true source. 

This directional-dependence of the angular flux proflle 
would be a smoking gun for pulsars causing the haze as 
opposed to the more spherical dark matter explanation. 
The Planck probe's increased sensitivity, larger number 
of bands, enhanced models, and expected smaller Galac- 




FIG. 4. Relative flux $ versus Galactic coordinates i and h for the standard NFW dark matter profile with 7 = I.O (left) and 
the NFW-type dark matter proflle with 7 = 1.2 (right). The mask is that from WMAP K-band mask Kp4 [4^ . 
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FIG. 5. Angular flux proflle, <1> versus angle from the Galactic 
center a for the 7 = 1.2 dark matter proflle and the Gaussian 
pulsar profile. The dashed line is the average over tl) — 30° — 
90° for the 7 — 1.2 dark matter distribution. The upper 
(dash-triple dotted) line and lower (dotted) line are the tp = 
30° and tp = 90° curves for the Gaussian pulsar distribution, 
respectively. The solid curve is an average from 7/) = 30° — 90° 
for the Gaussian pulsar distribution. The plots are scaled 
with the ip = 90° plot normalized to the flux of the haze. The 
points are the WMAP haze with errors from Ref. 0. 



tic mask will test these models [2SH27]. This will be in- 
strumental in determining the source of the WMAP haze 
as an astrophysical signal or the indirect detection of the 
dark matter. 
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